LOCAL EQUILIBRIUM OF THE QUARK-GLUON PLASMA 
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Within kinetic theory, we look for the local equilibrium configurations of a quark-gluon plasma by 
maximizing the local entropy. We use the well-established transport equations in the Vlasov limit, 
supplemented with the Waldmann-Snider collision terms. Two different classes of local equilibrium 
solutions are found. The first one corresponds to the configurations that comply with the so- 
called coUisional invariants. The second one is given by the distribution functions that cancel the 
collision terms, representing the most probable binary interactions with soft gluon exchange in the 
1/^ ' f-channel. The two sets of solutions agree with each other if we go beyond these dominant processes 

and take into account subleading quark-antiquark annihilation/creation and gluon number non- 
conserving processes. The local equilibrium state appears to be colorful, as the color charges are not 
' locally neutralized. Properties of such an equilibrium state are analyzed. In particular, the related 

^ . hydrodynamic equations of a colorful fiuid are derived. Possible neutralization processes are also 

' briefiy discussed. 
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^ ■ I. INTRODUCTION 

I ' In the course of equilibration a many-body system first reaches a local equilibrium and then it evolves hydrody- 
Oh, namically, usually at a much slower rate, towards global equilibrium. The distribution function of local equilibrium 
^ ■ is typically of the form of global equilibrium, but its parameters - temperature, hydrodynamic velocity, chemical 
potentials - are space-time dependent. However, the local equilibrium can also qualitatively differ from the global one. 
For example, the electron-ion plasma, which is homogeneously neutral in global equilibrium, can be locally charged 
before the global equilibrium is reached, see e.g. [1]. Thus, parameters that are irrelevant for global equilibrium might 
be needed to describe local equilibrium. While the state of global equilibrium is unique, the local equilibrium evolves 
' and even its qualitative features can change in time. The processes of charge neutralization are, for example, very fast 
in the electron-ion plasma. Therefore, the system is locally neutral after a short time but the electric currents survive 
for much longer. Thus, we deal with various local equilibrium states, depending on the time scale of interest. The 
form of local equilibrium is an important characteristics of a system. Knowing the respective distribution function, 
one can formulate a hydrodynamic description of the system. Let us again refer to the case of the electron-ion plasma. 
The fact mentioned above that the plasma is neutralized fast but the currents flow for a longer time justifies the 
magneto-hydrodynamics with no electric fields. 

The aim of this paper is to discuss local equilibrium of the quark-gluon plasma. While the global equilibrium features 
of the system have been studied in detail, see e.g. the review [2], not much is known about its local equilibrium. 
Although the problem was formulated long ago [3-6], the key questions remain unanswered. In particular, the 
scenario of equilibration of color degrees of freedom is far not established. It is unclear whether the regime analogous 
to magneto-hydrodynamics in the electron-ion plasma occurs in the quark-gluon plasma. However, the Yang-Mills 
magneto-hydrodynamics has been already considered [3,5-10]. 
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We intend to address these issues which are now of particular interest because of the large scale experimental program 
at the Relativistic Heavy-Ion Collider (RHIC) in Brookhaven National Laboratory, where high-energy nucleus-nucleus 
interactions are studied, sec e.g. [11]. At the early stage of such a collision, when the energy density is sufficiently high, 
the generation of the quark-gluon plasma is expected. The most spectacular experimental result obtained by now at 
RHIC is presumably an observation of a large magnitude of the so-called elliptic flow [12]. The phenomenon, which 
is just sensitive to the collision early stage, is naturally explained within hydrodynamics as a result of large density 
gradients [13]. Since the hydrodynamic description is applicable for a system in local thermodynamic equilibrium, 
the large elliptic flow suggests a surprisingly short, below 1 fm/c [14], equilibration time. Other characteristics of 
relativistic heavy-ion collisions are also consistent with a model assuming equilibrium state of strongly interacting 
matter produced in the collisions, see e.g. [15]. Thus, understanding of the equilibration mechanism of the quark-gluon 
plasma is a key problem for RHIC physics. 

The question of local equilibrium is related to a serious difficulty of the transport theory of the quark-gluon plasma. 
The local equilibrium is defined as a state which maximizes the local entropy. However, the entropy production 
occurs not due to the Vlasov evolution, which is rather well understood [2,16], but this is a dissipative phenomenon 
caused by the particle collisions. Thus, the collision terms of the transport equations are needed to discuss the local 
equilibrium. However, a derivation of these terms has occurred to be a very complex task and only the special case 
of quasi-equilibrium plasma has been seriously examined [17-21]. Fortunately, the structure of the collision terms 
can be guessed referring to the analogies between the spin and color systems. And this is not only a superficial 
similarity of degrees of freedom governed by the SU(2) and SU(3) group, respectively. The relationship appears to be 
much deeper. The covariancc of spin dynamics with respect to the rotation of quantization axis strongly resembles the 
gauge covariance of QCD. Thus, it was argued long ago [22] that the QCD collision terms are of the Waldmann-Snider 
type [23] known from the studies of spin systems. More recently, guided by the same analogy, the Waldmann-Snider 
transport equations have been used to compute color conductivity of the quark-gluon plasma [19], as well as other 
transport coefficients [24,25]. 

Once the collision terms of transport equations are known, the problem of finding the state of local equilibrium is 
well posed, see e.g. [26]. Namely, one looks for a configuration which maximizes the local entropy. In fact, such a 
configuration can be also found without a detailed knowledge on the collision terms. One only needs the so-called 
collisional invariants - the conditions obeyed by the coUision terms, coming from the conservations laws. In such 
an approach, already followed in [4,6], we, however, gain no information about the time scale corresponding to the 
local equilibrium sate. We also do not know whether the local equilibrium configuration dictated by the collisional 
invariants is the most general maximum entropy state. To answer these questions an explicit form of the collision 
terms is required. Then, one looks for a configuration that cancels the collision terms. 

In this paper we follow both approaches. After introducing the kinetic theory of the quark-gluon plasma in Sec. 
II, we find in Sec. HI the local equilibrium state provided by the collisional invariants. Then, we select the most 
probable binary interactions and we derive in Sec. IV the local equilibrium functions which cancel the Waldmann- 
Snider collision terms corresponding to these dominant processes. The derivation requires solving a whole set of rather 
complicated matrix equations. To simplify the analysis, we consider particles obeying classical statistics, although 
we believe that the physical picture emerging from our analysis is not much changed when quantum statistics is 
incorporated. The local equilibrium states, which come from the approaches of Sees. HI and IV, are colorful and their 
color structure is exactly the same. However, the baryon chemical potential of (anti-)quarks and the scalar chemical 
potential of gluons remain unconstrained by the dominant processes. The constraints provided by the collisional 
invariants only appear when the subleading quark- antiquark annihilation/creation and gluon number non-conserving 
processes are included. To better understand properties of the colorful local equilibrium, we derive in Sec. V the 
resulting hydrodynamic equations. Finally, we consider the applicability of our results and briefly discuss possible 
processes responsible for the color neutralization in the quark-gluon plasma. Some formulas of the S\]{Nc) generators 
are collected in the Appendix. 

Throughout the paper (except Eqs. (115-117) where c is restored) we use the natural units with c = U = ks = ^ 
and the metric (1, —1, —1, —1). 

II. KINETIC THEORY OF THE QUARK-GLUON PLASMA 

In this section we discuss the transport theory of quarks and gluons [16,27]. The SU(Afc) gauge group is left 
unspecified but we pay a particular attention to the cases Nc = 2 and Nc = 3, for their possible applications to the 
different high temperature phases of the Standard Model. Generically speaking, we call gluons the particles associated 
to the vector bosons of SU(A^c)j which carry charge in the adjoint representation, and we call quarks or antiquarks 
the particles with the charge in the fundamental representation. We will also call parton any of those particles. 
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A. Distribution functions and transport equations 



The distribution function of quarks Q{p, x) is a hermitian x matrix in color space (for a SU(A'^c) color group); 
X denotes the space-time quark coordinate and p its momentum, which is not constrained by the mass-shell condition. 
The spin of quarks and gluons is taken into account as an internal degree of freedom. The distribution function 
transforms under a local gauge transformation U as 



Q{p,x) U{x)Q{p,x) U\x) 



(1) 



that is, it transforms covariantly in the fundamental representation. Here and in the most cases below, the color 
indices are suppressed. The distribution function of antiquarks, which we denote by Q{p,x), is also a hermitian 
Nc X Nc matrix in color space, which in a natural way should transform covariantly in the conjugate fundamental 
representation. However, we will express the antiquark distribution function in the same representation as quarks 
throughout, and then it transforms according to Eq. (1). The distribution function of (hard) gluons is a hermitian 
(A^^ — 1) X (A^^ — 1) matrix, which transforms as 



where 



with ■ 



G{p,x)^l{{x) G{p,x) U\x) 



Uab{x) = 2Tr[T"t/(x)T''t/+(x)] , 
, — 1 being the SU(A^c) group generators in the fundamental representation with Tr(r'*T'') — ^ 



We note that =U ^ = W . Therefore, not only G but also G^ transforms covariantly i.e. 

G^{p,x) ^U{x) G'^{p,x)U\x) . 
The color current is expressed in the fundamental representation as 



where the momentum measure 



■1/ 



dP p>' [Q{p, x) - Q{p, x) - ^Tr[Q(p, x) - Q{p, x)] + 2t«TV [T«C;(p, x)] 



dP 



d'^p 
(2^ 



2Q{po) 5(p2 



(2) 

(3) 

(4) 
(5) 

(6) 



takes into account the mass-shell condition po = |p|. Throughout the paper, we neglect the quark masses, although 
those might be easily taken into account by modifying the mass-shell constraint in the momentum measure. A sum 
over hclicitics, two per particle, and over quark flavors Nf is understood in Eq. (5), even though it is not explicitly 
written down. The SU(iVc) generators in the adjoint representation are expressed through the structure constants 
Tj" = —ifabc- and are normalized as Tr[T"T^] = NcS"'''. The current can be decomposed as j'^{x) = j^{x)T'^ with 

j^{x)=2TT{TaJ^{x)). 

Gauge invariant quantities are given by the traces of the distribution functions. Thus, the baryon current and the 
energy-momentum tensor read 



h^'{x) = ^jdPp''Ti 



Q{p,x) - Q(p, x) 
''{x)= j dP pf'p" Tr [q(p, x) + Qip, x) + G{p, x) 



(7) 
(8) 



where we use the same symbol Tr[- • •] for the trace in the fundamental and adjoint representations. 
The entropy flow is defined as [4] 



s^'{x) = - j dPp^'Tr 



QlnQ + (1 - Q)ln(l - Q) + QlnQ + (1 - Q)ln(l - Q) 

+ G\nG -{1 + G)ln(l + G) 



(9) 
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If the effects of quantum statistics are neglected, Eq. (9) simplifies to 

s^'{x)=- j dPp^'Tr Q{\nQ-l) + Q{lnQ-l)+G{\nG-l) . (10) 
The distribution functions of quarks and gluons satisfy the transport equations: 

j/^D^Q{p,x) + ^j/'{F^.{x),d;Q{p,x)} =C[Q,Q,G], (11a) 

p^^D^Q{p,x) - ^p^- {F^,{x),d;Q{p,x)} =C[Q,Q,G], (lib) 

P^V^Gip, x) + ^p^ {J'^Ax), d;G{p, x)} = Cg[Q, Q, G] , (11c) 

where g is the QCD coupling constant, {..., ...} denotes the anticommutator and dp the four-momentum derivative; 
the covariant derivatives and "D^ act as 

D^=di^- ig[A^{x), ... ] , ^^=81^- ig[A^{x), ... ] , 

with An and A/j. being four-potentials in the fundamental and adjoint representations, respectively: 

A'^ix) = ^^(a;)T« , ^"(a;) = TMJ^(x) . 

The stress tensor in the fundamental representation is F^i, = d^A^, — d„Afj_ — ig[Afj_, A^], while T^,, denotes the field 
strength tensor in the adjoint representation. The collision terms G,C and Gg are discussed in detail in the next 
subsections. 

Let us finally mention that in the transport theory framework one can consider two different physical situations: 

1) the gauge fields entering into the transport equations (11) are external, not due to the plasma constituents; 2) the 
gauge fields can be generated self-consistently by the quarks and gluons. In the latter case, one also has to solve the 
Yang- Mills equation 

DnF^''{x)=f{x), (12) 

where the color current is given by Eq. (5). 



B. Decomposition of the distribution functions and associated transport equations 

The parton distribution function N is essentially the statistical average of the Wigner transform of the product 
of two field operators representing quarks or gluons [16]. If the parton carries color charge in a representation R, 
then the distribution function N transforms under gauge transformations as R, where R is the representation 
conjugate to R. 

In the SU(2) group, the products of the fundamental (2) and adjoint (3) representations decompose into irreducible 
representations as 

2®2==le3, (13) 
3(g)3 = l©3©5. (14) 

As known, the conjugate and direct fundamental representations of SU(2) are equivalent to each other. The decom- 
position of the products of the fundamental (3) and adjoint (8) representations of the SU(3) group are 

3®3 = 1®8, (15) 

8®8 = ie8e8eiOBToe27. (i6) 

The above decompositions show that the distribution functions of quarks and antiquarks are uniquely specified by 
their singlet and adjoint components. Thus, the functions can be written as 

Q{p,x) = ■^qo{p,x)+q''{p,x)T'' , (17a) 

Q{p, x) = -^qoip, x) + (t{p, x)t'' , (17b) 
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where 

qo{p, x) = Tr[Q(p, x)] , qa{p, x) = 2TrKQ(p, x)] , (18a) 
qo{p, x) = Tr[Q{p, x)] , q^ip, x) = 2Tr[T'^g(p, x)] . (18b) 

From Eq. (11a) it is possible to deduce a set of coupled equations for the colored and colorless components of the 
quark distribution function which read 



P^d^p, x) + I p^F^^ix)^^!^^ = Tr[C] , (19a) 



p'^D;\{p,x) + I d.,,p^F^^^i^)^^^ + j_p.F;^{x)^-S^ = 2Tr[r'^C] , (19b) 

where date are the totally symmetric structure constants of SU(A^c) and D'^'^ = 9^(5'"^ + gfatcA''^- The projected equa- 
tions, which can be also written for antiquarks, show that transport phenomena of colorless and colored components 
are coupled beyond the lowest order in the gauge coupling constant. 

From the decompositions (14,16) it is clear that the singlet and adjoint components arc not enough to fully de- 
scribe the gluon distribution function. For gluons one also needs components in higher dimensional representations. 
Below, we present a way to uniquely characterize the gluon distribution function in terms of its fully symmetric and 
antisymmetric components for the SU(2) gauge theory. 

We first express G{p, x) as 

G{p,x)^gab{p,x)T''T\ (20) 

which uses as a basis for 3 x 3 hermitian matrices the set of 9 independent matrices T°'T^. We note that both G and 
Q are 3x3 matrices which are related to each other as 

Gab{p, X) = d'^^'GUP, X) - GbaiP. ■ (21) 

Expressing the product of T"T'' as 

yayb = i [T", T^] + ^{T", r''} , (22) 
and taking into account that the commutator is proportional to T'^, instead of Eq. (20) we write 

G{P, x) = ^gaip, x) T'^ + gabip, x) ^{T% T'} , (23) 



where 



ga{p, X) = iP^" Qcbip, X) , gab{p, x) = ^ {Gab{P, X) + QbaiP, x)) . (24) 



The equation (23) can be also written as 

GabiP, x) = -ir'"'gc(p, x) + 5"^5cc(p, x) - gabip, x) . (25) 

Thus, according to the decomposition in Eq. (14), the antisymmetric components of G correspond to the representation 

3, while the 6 symmetric components correspond to the 5 and 1, the last one being the trace. Because of the Casimir 
constraint, TT" = 2, the singlet component can be obtained from the symmetric part gab 

go{p, x) = Tr[G(p, x)] = 2 gaa{p, x) . (26) 

The transport equations obeyed by Qq, ga and gab are found multiplying Eq. (11c) by the unity, T° and {T°-,T'^}/2, 
respectively, and taking the trace. Using the relations (A12), we get 

p'^dMP^ x) + gpi^F^^ix)^^^^ = Tt[C,] , (27a) 
p^Dfgbip, X) + gp^Fl,Ax) + ^^^^) = ^^[T'^C,] , (27b) 

P^{D,)llg^M + ^jf (^J.;^(x)^^^ +F^,(x)^^^) = \tv[{T^T'}C,\ - \6^'T,[C,] , (27c) 
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where 

{D^)t^ = d^S'^'^S^'' + gn^S'"' Al + gf^H'^'^ (28) 

is the covariant derivative acting on a tensor of rank 2. Note that multiplying the last equation by ^""^ ^ we get, as 
expected, the equation for go- 

For SU(3), or SU(A^c) in general, the decomposition of the gluon distribution function into irreducible representations 
and the equations obeyed by every component have a much more involved structure, and they will not be discussed 
here. 



C. Waldmann-Snider Collision terms 



The transport equations for the quark-gluon plasma (11) have been written down without specifying the collision 
terms. Unfortunately, a complete derivation of C, C and Cg is still lacking, as already mentioned in the Introduction. 
However, using the analogy with the spin systems one can justify the use of the Waldmann-Snider collision terms. 
The main characteristic of these collision terms is that they depend on the scattering amplitudes rather than on the 
collisional cross sections, as it happens in the usual Boltzmann equation. 

Let us discuss the general structure of a collision term for a system of particles carrying quantum color charges. 
The most probable processes are binary collisions (p, r;pi.s) ^ {p' ,f/,p'i,u) where p,pi,p' ,p'i denote the momenta 
and r, s, t, u colors, in the fundamental or adjoint representation, of interacting partons. We denote by N{p, x) the 
generic distribution function of the partons - quarks or gluons. The Waldmann-Snider collision term, which enters 
the kinetic equation of N, is of the form [26]: 



C[N, Nu N', Ni] = J dP'dPidP, (27r)^^(^) (p + p,-p'- p[) [i {l ± N, J+} - 1 {N, !_} 



(29) 



where we have used a rather common notation N = N{p,x), Ni = N{pi,x), N' = N{p',x), and = N{p[,x). The 
first term, which represents a gain term, is given by 

J7 = Mrstu{p,Pi;p',p[) M;,^^{p,p,;p',p[) N'\p',x) N^Hp[,x) {l±N{p^,x)r , (30) 

while the second one is a loss term defined as 

IL' = Mrstu{p,Pi;p',p[) M;,^^{p,p,-p',p[) 7V-(pi,x) (1 ± N{p',x)f (1 ± N{p[,x)r . (31) 

Mrsiu represents the scattering amplitude associated with the collision process under consideration. The double sign 
± reflects the fermionic character of quarks and bosonic of gluons. We have used here the compact notation of ref. 

For the consistency of the theory, it is necessary to prove that the Waldmann-Snider collision terms transform 
covariantly under a gauge transformation, in the same way as the left hand sides of the transport equations (11) 
do. It is difficult to check this gauge covariance in full generality without specifying the scattering process and 
the corresponding scattering amplitudes. For all the cases we are going to consider, the gauge covariance of the 
Waldmann-Snider collision term holds as the distribution functions transform covariantly (see Eqs. (1) and (2)), and 
the scattering amplitudes, stripped of the color generators, are gauge invariant. We will briefly come back to this 
point in Sec. IV. 



D. Conservation laws and entropy production 



As well known, the collision terms should satisfy certain relations due to the conservation laws. In our case, the 
laws are: the baryon charge conservation 

d^bi^ix) = , (32) 

the energy-momentum conservation 

d^t^'-'ix) + 2Tr[j,{x) F^-'ix)] = , (33) 
and the covariant conservation of the color current 
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D^j^ix) = . (34) 

Let us derive the relations constraining the coUision terms which follow from Eqs. (32,33,34). Using the transport 
equation, one finds from the definition (7) 

d^b^ix) = lJdPTr[C-C]-^J dPp'^ Tr[F^.d;iQ + Q)] . 

Now, one performs partial integration of the second term in the r.h.s. Assuming that the distribution functions vanish 
at infinite momentum and observing that g'^^Ffj^i, = 0, one finds that the term equals zero. Therefore, the baryon 
current conservation (32) provides 



dPTr[C-C] =0. (35) 
In analogous way, one finds that the energy-momentum conservation (33) implies 

dP Tr[C + C + Cg]=0, (36) 
while the covariant conservation of the color current leads to 

dP[C-C + 2r"Tr[r°Cg]] = , (37) 



where we have taken into account the relation (35). 

Let us now discuss the entropy production. We neglect here the efi'ects of quantum statistics, and consequently 
start with the definition (10). Following the derivation of Eqs. (35,36,37), one finds 

di^s^ix) = - j dP Tr[ClnQ + ClnQ + CglnG] (38) 

dP p^ Tr [{F^,, Q}d;\nQ - {F^,, Q}d;\nQ + {T^,, G}9;inG] , 
where the partial integration has been once performed and it has been observed that 



Tr 



= 



and that the analogous equalities hold for Q and G. Assuming that Q and d^Q commute with each other i.e. 

[Q,5;O]=0, (39) 

one shows that d^hiQ = Q~^dpQ. Using the condition (39) and the similar ones for Q and G, one proves that the 
second term in r.h.s. of Eq. (38) vanishes after one more partial integration. Then, we get 

d^s^'ix) = - j dP TrfClnQ + ClnQ + CglnG] . (40) 

According to Eq. (40), the entropy of the quark-gluon system is produced due to the collisions. If the commutation 
condition (39) is relaxed, the second term in r.h.s. of Eq. (38) does not vanish, and we arrive to a paradoxical result 

that the mean-field dynamics docs not conserve the entropy. 

A local equilibrium configuration is achieved when there is no entropy production, i.e. 9^s^(x) = 0. This equation 
is of very complicated structure and it has two classes of solutions. The first one cancels the collision terms but to 
get it the collision terms have to be specified. The second class appears due to the conservation laws, i.e., because of 
the relations (35,36) and (37). In the remaining part of this article, we will study the two sets of solutions. 
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III. LOCAL EQUILIBRIUM FROM THE CONSERVATION LAWS 



In this section we discuss, following [4,6], consequences of the conservation laws (32,33,34). Specifically, we obtain 
the local equilibrium configuration which is found as a solution of the equation 



/ 



dP Tr[ClnQ + ClnQ + CglnG] = 



(41) 



due to the relations (35,36,37). 

One easily constructs the local equilibrium distribution function out of the collision invariants. Indeed, one shows 
using Eqs. (32,33,34) that Eq. (41) is solved if 



Qeq{p,x) = exp - /3{x){u„{x)p'' - iib{x) - /I(x)) 
Qeq{p,x) = exp - l3{x){u„{x)p'' + IJ.b{x) + Jiix)) 
Geq(p, x) = exp - P{x){u„{x)p'" - Jlg{x)) 



(42a) 
(42b) 
(42c) 



where /3(.x), u'^(x) and fibix) arc, respectively, the inverse temperature, hydrodynamic velocity and baryon chemical 
potential which are all scalars in color space. The color chemical potentials Ji and Jig are hermitian matrices Nc x Nc 
for quarks and {N^ ~ 1) x (-^c ~ 1) ^'^^ gluons. They are gauge dependent variables, which transform as 



Jlix) ^ U{x) Jl{x) uHx) , /7g(;r) ^ U{x) /7g(.T) U^x) 



(43) 



In general, Ji can be expressed as Jl = fxo + HaT"'- However, the singlet component /xq is already singled out as a 
baryon chemical potential /Ufo. Therefore, we write down ju = /Uar". Consequently, the color chemical potential /I is 
not only hermitian but also traceless. The covariant conservation of the color current provides the relation 



Jig = 2T»TV[t"m] = MaT" , 



(44) 



which implies that Jig is also traceless. The baryon and color chemical potentials occur in Eqs. (42) because of the 
conservation laws of baryon number and color charge, respectively. The temperature and hydrodynamic velocity are 
related to the energy-momentum conservation. 

The local equilibrium state described by Eqs. (42) is not color neutral. Substituting the distribution functions (42) 
into Eq. (5) one finds the color current as 



^3 



^/(e^'"'(e'^'^- lJ-Tr[e^'']) 



+ 2T"Tr[r"e'^^3] 



(45) 



where T, u^, /i^, Jl, and Jig are functions of x. The fact that the color current is finite does not imply that the system 
as a whole carries a finite color charge. We note that the a;— dependence of the color chemical potentials, which enter 
the solutions (42), is not specified. Therefore, it can be always chosen in such a way that the total color charge defined 

as / cPxj^ vanishes. 

The derivation of the local distribution function based on the collisional invariants tells nothing about the time 
scales when the colorful configuration (42) exists. To get such an information the collision terms have to be specified. 
This is discussed in the next sections. 

The equilibrium solutions (42) are given in an arbitrary gauge. It is often useful to work in a gauge where the quark 
and antiquark chemical potentials are diagonal. Then, 



11 = H T 



drpd 



Jig = li^T 



(46) 



where r"^ and T'^ are the fundamental and adjoint generators of the Cartan subalgebra of SU(Afc) (d = 3 for SU(2) 
and (i = 3, 8 for SU(3)). In this gauge one has, as will be seen below, well-defined numbers of quarks and antiquarks 
of a certain color. And then, the physical meaning of the color chemical potentials becomes transparent. 
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A. Diagonal gauge for SU(2) 



Using the explicit form of = where is the Pauli matrix, the singlet and the (non-vanishing) adjoint 

components (see Eqs. (17,18)) of the quark and antiquark distribution functions of local equilibrium (42) are found 

as 



qo{p, x) = (71- (p, x) + qi (p, x) , 

qo {p, x) = (p, x) + qi (p, x) , 
where the scalar functions q^i and q^i are 



93 (P, x) = (p, x) - qi (p, x) , 
93 {P, x) = 9T (P> x) - qi (p, x) , 



-I3{x) ( u^{x)pij, - Ubix) T 2 '"3(2;) 



9Tib.a^) = exp 



-P{x) ( u''{x)pi^ + iJ,b{x) ± -M3(a;) 



(47a) 
(47b) 

(47c) 
(47d) 



While the generator is diagonal, is not. To derive the expressions for gluons one has to observe that [T^Y is 
the diagonal matrix with 1, 1,0 on the diagonal. Consequently, (T^)" = when n = 1,3,5 .. . and (T^)" = {T^Y 
when n = 2, 4, 6 Thus, the non- vanishing components (23) of the gluon distribution function (42c) are 

go{p, x) = g^{p, x) + g^{p, x) + ^^(p, x) , (48a) 



53 (p, x) = g^{p, x) - f;^(p, x) , 



511(^,3;) =g22ip,x) = -g^{p,x) 



533 (p, x) = ^ (a-ttip, x) + gi^ip, x) 

where the functions g^jj, and are 

mw^ip, •^) = cxp [-li{x) {u''{x)pf, =F ^3(2:))] , 



^(P,x)j , 



g^(p,x) = cxp [-/?(.T)?i^(x)p^ 



(48b) 
(48c) 

(48d) 

(48e) 



In the diagonal gauge, a finite value of the color chemical potential simply means that the populations of quarks, 
antiquarks and gluons of difl:crcnt colors are not the same. 

Using the distribution functions in the form (47,48), the color current (45) can be written as 



2~'3 



NfchiPut) sh(/JM3/2) + sh(/3At3) 



(49) 



B. Diagonal gauge for SU(3) 

The local equilibrium solutions for the SU(3) plasma can be also written in the diagonal gauge. However, the 
formulas are not that simple as for the SU(2) case. We take the generators in the fundamental representation as 
r" = A"/2, where A° are the Gell-Mann matrices. The matrices and A^ are diagonal with the elements 1, —1, 
and l/-\/3, 1/v^, —2/y/3, respectively, along the diagonal. With a color chemical potential in the directions a = 3 
and a = 8 one can then easily evaluate the singlet and (non-vanishing) adjoint components of Qeq and Qeq, which we 
write in terms of the distributions functions of red, blue and green quarks and antiquarks. Here, we have taken the 
convention to assign the first, second and third rows/columns of the Gell-Mann matrices to the red, blue and green 
colors, respectively. A simple evaluation leads to 

qo{p, X) = gred(P, X) + qUueiP, X) + 9green(p, x) , qi{p, x) = 9red(p, x) - qb\ue{p, x) , (50a) 

qs{p, x) = (^q^ed{p, X) + qUueiP, x) - 2ggreen(p, xfj , (50b) 

go(p, x) = qrediP, x) + qUue{P, x) + ggreen(p, x) , Q3(p, x) = qred{P, x) - qh\ue{p, x) , (50c) 

QsiP, ^) = [qrediP, X) + qh\ue{P, x) - 2qgreen{p, xfj (50d) 
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where the distribution functions of quarks and antiquarks of different colors are of the form (47c) and (47d), respec- 
tively, but with the following color chemical potentials: 

Hred{x) = - y^six) + j , tJ.Uue{x) = --yisix) - j , ^^green{x) = (51) 

The computation of the singlet and adjoint components of the local equilibrium distribution function of gluons 
is much more involved. The evaluation of the traces requires to expand the exponentials, and compute the traces 
of arbitrary powers of T3, of Tg, and of TaTg. With the help of Mathematica, we have found the singlet and (non- 
vanishing) adjoint components as 

goip.x) = 2g^{p,x) +g^+{p,x) +g^^(p,x) +gy+{p,x) +gy-{p,x) +g^+{p,x) +g^-{p,x) , (52) 

53 ip,x) = 3z+ {p,x) - g^-{p,x) + - (^g^+ (p, x) + g^-{p,x) + gy+ {p,x) + gy-{p,x)^ , (53) 

58 (P, ^) = -Y (^'^+ ^) ~ 9^- (P' ^) - 9y+iP,x) + gy-{p,x)^ , (54) 

where the scalar functions gs, 5x±, 9y±, and gz± are analogous to those from Eqs. (48e) but their color chemical 
potentials are 

Hs{x) = 0, Mx±(a;) = ±— ^± , /xy±(a;) = ±— ^ T , i^^±{x) = ±ii3{x) . (55) 

Exactly as in the SU(2) case, we find that a finite value of the color chemical potential means that quarks, antiquarks 
and gluons of different colors have different densities. 

IV. LOCAL EQUILIBRIUM FROM VANISHING COLLISION TERMS 

As follows from Eq. (40), there is no entropy production when the collision terms vanish. Thus, local equilibrium is 

reached when the gain and loss terms compensate each other. Consequently, we will look for solutions of the equation 
C — 0. However, there are numerous scattering processes occurring in the quark-gluon plasma and. in general, the 
complete set of collision terms entering into the quark, antiquark and gluon kinetic equations is rather large, even so 
we only consider the binary collisions. The most probable processes, i.e. those with the largest cross section, occur 
when two partons exchange a soft gluon in the t- or u-channels. The later possibility only happens for interaction 
of identical partons - quarks of the same flavour or gluons. In vacuum, the corresponding cross sections diverge as 
or when the four-momentum transfer t ov u goes to zero. In the medium, these divergences are softened, as 
the gluon propagators are dressed by the interactions, and the electric and magnetic forces are either statically or 
dynamically screened. In the local equilibrimn state, which is achieved at the shortest time scale, the collision terms 
associated with those processes, we call them 'dominant', have to vanish. Thus, we will first consider the interactions: 
qq qq, qq <-> qq, qq qq, gg <-> gg, qg <-> qg, and qg qg, and we will neglect all other processes, as they are 
relevant for longer time scales"'^. These less probable processes drive the system either to a different local equilibrium, 
or to the global equilibrium. We will also consider the subdominant processes with the soft quark in t- or w-channel 
which correspond to the vacuum cross sections diverging as or respectively. These are the quark-antiquark 
annihilation and creation into and from two gluons in t- or u-channel which, as will be shown, have a qualitative 
effect on the local equilibrium state. With the subdominant processes, one should also consider all the channels and 
the respective crossing tcirms of the various binary collisions, plus another set of collisions that do not conserve the 
particle number. The complete analysis is very complex, and we will not carry it out here. 

In this section we write down the relevant collision terms, and then we discuss the equations imposed by the 
vanishing of these terms. Finally, we solve the equations, showing that the nature of local equilibrium is fixed by the 
color structure of the scattering amplitudes. 



^To estimate a mean free time associated with a given collision process one has to specify the distribution function. For a 
discussion of those mean free times in global equilibrium see Ref. [28] . 
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A. Collision terms 

The dominant parton-parton scattering amplitudes with one-gluon exchange in the t- and u-channels are of the 
form 

Mrsr's'{p,Pi;p',p'i) = M{p,pi;p',p[) T^»,7;«' (56) 

Mrsr's'{p,Pi;p',p[) = M{p,Pi;p\p[) Vs'Vr' , (57) 

where T" and T" arc the group generators of SU(A^c) of the two partons participating in the colhsion: T" — T° 
for gluons, T" ~ t" for quarks and T° — — (r")"^, where T means transposition, for antiquarks. With the t-channel 
amphtude (56), the colhsion term (29) equals 

C[N, iVi, N', N[] = J dP'dP[dPi (27r)^5W ip + Pi-p'- p'l) \Mf (58) 

X (t^N'T^ Tr[f''N[f^] - ^{T^T",7V} Trft^TVit'']) , 

where we have neglected the effects of quantum statistics, and consequently the terms 1± N have been replaced by 
unity. The collision term corresponding to the 7i-channel amplitude (57) can be found from (58) by means of the 
exchange N -i^ Ni and N' ^ in the r.h.s of Eq. (58). 

Using the identity (Al) given in the Appendix, we can write down Eq. (58) for the case of quark-quark scattering 

as 

C[Q,Qi,Q',Q[] = dP'dP[dPi (27r)V4)(p + pi -p'-p'^) \M\^ (59) 

X ((Tr[Q']g; - Tr[Q]gi) - -i^(Q'Tr[Qi] - QTrfQi]) - ^{{Q'M - {QM)) • 
The collision term (58) for the quark-antiquark scattering is 

C[Q,Qi,Q' ,Q[] = dP'dP[dPi (27r)^(5(^)(p + pi - p'l) \M\'' (60) 

X ((Tr[Q'(5'i] - ^{Q,Oi}) - -^(0'Tr[Q;] - QTr[gi]) - j^{{Q' ,Q[} - {Q,Qi})) , 

where, as discussed previously, we have replaced by Q. 

For the gluon-gluon scattering we have found a simplification of Eq. (58) only in the case of the SU(2) gauge group. 
Then, the collision term reads 

C[G,Gi,G',G"i] = j dP'dP[dPr {2n)^d^^\p + - p' - p[) {Mf (61) 

X ((Tr[G'^G;] - {G'^,G[} - i{G,Gf}) + (Tr[G'] G[ - G Tr[Gi])) . 

The scattering amplitudes of the subdominant processes with the quark exchange in t- and u-channel have the 
following color structure 

M^jab{p,Pi;p',p'i) = M{p,pi;p' ,p\Kk4j , (62) 
Mijab{p,Pi\p' ,p'i) = M{p,pi;p',p[)TfkT^j . (63) 

The collision term associated with this f-channel annihilation processes is 

G[g,Qi,G',G'i] = J dP'dPidPi (27r)*<5W(p + pi -p' -p'l) \Mf (64) 

X (tVVMg««(p')G''"^(p'i)-^{Q(p),tV^Q(Pi)tV«}) . 
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At the end of this section we call the attention of the reader to the structure of the collision terms (59,60, 61). 
Because there are only objects like Q, G, , which transform covariantly with respect to the gauge transformation 
(1,2), and Tr[(5], Tr[G] and Tr[Q(5i], which are gauge invariant, these collision terms transform covariantly, provided 
\M.\^ is gauge invariant. The gauge covariance of the collision term (64) is evident when instead of C the projections 
Tr[C] and Tr[T'*C] are considered. As will be seen in the following subsection, these projections have the right gauge 
structure. 



B. Conditions of local equilibrium 

In this subsection we present the conditions for the cancellation of the collision terms associated with the processes 
discussed above. 



1. qq <-> qq 

The collision term (59) corresponding to the quark-quark scattering vanishes if 

(Tr[0']Q; - Tr[g]Oi) - 1^{Q'Tv[Q'^\ - QTv[Q^]) - ^{{Q' ,Q{} - {Q,Qi]) = , (65) 

where p + pi =p' +p'\- Because the quark matrix transport equation can be uniquely characterized by its singlet and 
adjoint components (see Eqs. (19)), the condition (65) requires 

Tr[QQi] = Tr[Q'g;] , (66a) 

TV[Q] Tr[gi] = Tr[Q'] Tr[Q;] , (66b) 

and 

Tr[T« {Q, Qi}] = Trir'^ {Q', Q'J] , (66c) 

Tr[r»Q] Tr[Qi] = Tr[Q'] Trfr^Qi] , (66d) 

TV[r»Q] TV[Qi] = Tr[r»Q'] Tr[Q;] . (66e) 

The conditions for cancellation of the collision term for antiquark-antiquark scattering are totally analogous to 
those of the quark-quark case. 



2. qq <-> qq 

The collision term (58) for the quark- antiquark scattering vanishes when 

(Tr[Q'g'i] - ^{g,gi}) - ^{Q'T:^[Q'A - QTr[gi]) - 1-{{Q',Q'^} - {Q,Qi}) = . (67) 

The conditions of cancellation of the projected matrix equation (67) read 

Tr[Q(5i] = Tr[Q'Q[] , (68a) 

Tr[Q] Ti[Qi] = Tr[Q'] Tr[Q[] , (68b) 

and 

Tr[T'^{g,Qi}] - Tr[T«{g',gl}] = 0, (68c) 

Tr[T"g] Tr[gi] = Tr[T"g'] Tr[g'i] . (68d) 

The requirement that Tr[T'' {g, gi}] = directly follows from the first term of Eq. (67). 
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3. qg <-> qg 



The collision term for quark-gluon scattering with one-gluon exchange in the t-channel vanishes if 

r'^gV' TrlT^G'iT''] - ^{r^'r", Q} Tr[T"Gir''] = . (69) 

Requiring that Tr[C] =0 and Tr[T"C] = provides the equations 

Tr[Q] Tr[Gi] - ^Tr[r'^Q] Tr[T"Gi] + irf"''^ Tr[r"Q] Tr[T*'T=Gi] 
= Tr[Q'] Tr[G'] - ^Trfr^Q'] Tr[T°Gi] + \cP^'= Trfr^Q'] Tr[T^T^Gi] , (70a) 

and 

Tr[T'=T''QT"] Tr[r''Gir''] + R'lQ, Gi] = T^It^t^Q't"-] T^lT^G'^T"] , (70b) 



where 



Gi] = |/cad(TV[T'^T*Q] TV[T*GiT«] - Tr[rV''Q] Tr[T«GiT*]) . (70c) 



4- 99 ^ 99 



The collision term of the gluon-gluon scattering equals zero when 

T^G'T^ Tr[T"G'iT''] - ^{T^T'^, G} Tr[r''GiT''] = . (71) 

For the SU(2) plasma the above condition can be simplified (see Eq. (61)) and it gives 

(Tr[G'^G;] - {G'^, G'J - \{G, Gf }) + (Tr[G'] G'^ - G Tr[Gi]) = . (72) 

We demand the cancellation of the totally symmetric and antisymmetric components of (72), see Eqs.(27). Imposing 
Tr[TaGg] = and Tr[{ra, TfcjGg] = 0, we get 

Tr[Ta {G, Gf}] = Tr[r„ {G', Gf }] = , (73a) 

Tr[r„G] Tr[Gi] = Tr[G'] Tr[T,G;] , (73b) 

and 

TV[{T„, n} {G, Gl]] = 85'^'' Tr[G'Gf ] - 2Tr[{T„, n} {G', Gf }] (73c) 

Tr[{T„, n}G] Tr[Gi] = TV[G'] Tr[{T„, Tb}G[] . (73d) 

For b Eq. (73c) requires that 

Tr[{T„, n} {G, Gj}] = Tr[{Ta, U} {G', Gf }] = , (73e) 

while for a = 6 (r"T" = 2) Eqs. (73c,73d) imply 

Tr[GG^] = Tr[G'Gf ] , (73f) 
Tr[G] Tr[Gi] = Tr[G'] Tr[G;] . (73g) 
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5. qq^gg 



With the scattering ampHtude given in Eq. (62), the cancellation of the collision term corresponding to the quark- 
antiquark annihilation in the t-channel demands 

T^It'^Qt'^t^Qxt^] = Trlr^T^T^T^] G"'^G'f , (74a) 
^Tr[r« {Q, t^t^'QitV"}] = Tr[rV'=r''TV] G'^^G'^ . (74b) 

The left hand side of the above equations can be simplified using the relations (A4) given in the Appendix and the 
formula 

1 _ _ 9 

t'^t'Q^t^'^ = + -^T^m , (75) 

Furthermore, for the SU(2) plasma one finds, using the relations (A7, A8) given in the Appendix, that 

Tr[QQi] + 4Tr[g] Tr[Qi] = 2Ti[G'^G[] - 2Tt[G'G[] + 2Tr[G'] Tr[G'i] , (76a) 
Tr[T^ {Q, Qi}] + 8Ty:[T^Q] Tr[Qi] = 2G'^«Tr[T'=G;] + 2Tr[T^Gf G'] - 2Tr[T^G;G'] + 2Tr[T'=G'] Tr[G;] . (76b) 



C. Local Equilibrium Solution for the SU(2) plasma 

Wc find here the local equilibrium solutions that cancel all the collision terms discussed in the previous subsection 
for the SU(2) plasma. We start with the quark-quark scattering. Eqs. (66a,66c) are solved by functions obeying 

Q{p, x) Q{pu x) = Q{p', x) Q{p[, x) , (77) 

for p + pi = p' + p'l- Using standard arguments, see e.g. [26], one finds that Eq. (77) is satisfied by exponential 
functions 

Q{p, x) = exp - I3{x){u^{x)p'' - iJ.b{x) - Ji{x)) , (78) 

where u^{x) and fl{x) are hermitian and traceless matrices. Please note that the scalar chemical potential Hb, which 
is interpreted as the baryon chemical potential, is already singled out. Because of Eq. (39), u^{x) and fl{x) should 
obey the condition [m^(x), /i(x)] — 0. Thus, using the gauge freedom to rotate these quantities in color space, they 
can be chosen in diagonal form. 

Eqs. (66b,66d,66e) require that the hydrodynamic velocity u'^{x) is proportional to the unit matrix. Otherwise 
different components of u^{x) enter differently Eqs. (66b,66d,66e) and the constraint p + pi — p' + p'l is insufficient 
to satisfy these equations. Once u'^{x) is proportional to the unit matrix, the condition ['u^(.t), /i(a;)] = is trivially 
satisfied, and there is no reason to require jl{x) to be diagonal. It is then an arbitrary traceless matrix, even so it 
can be diagonalized because of the gauge freedom. Since the hydrodynamic velocity is no longer a color matrix but 
a scalar, it is from now on denoted as m'' not as 77/*. 

Repeating fully analogous considerations for the collision term of antiquark-antiquark scattering, we arrive to the 
antiquark distribution function 

Q{p,x)=exp - (3{x){u^{x)p'' + j2b{x) + il{x)) . (79) 

The conditions of cancellation for the quark-antiquark collision term provide the relations between the parameters 
of quark and antiquark distribution functions. Namely, Eqs. (68a,68b,68d) require 

f3{x) u^ix) = P{x) u^ix) . (80) 

Because u''(a;)u^(x) = u''(a;)u^(x) = 1, we effectively have 

u^{x) = u^'ix) , T{x) = f{x) . (81) 

Furthermore, Eq. (68c) imposes 
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fl{x) = -Ji{x) , (82) 

but it leaves the baryon chemical potentials fXb and p,b unrestricted. 

Let us find now the distribution functions that cancel the gluon-gluon collision term. Condition (73f) is solved by 
those functions obeying 



which demands that 



G{p,x)G^{pi,x) = G{p',x)G^{p[,x) , (83) 



G{p,x) = exp - Pg{x){u^g{x)p^ - ijPg{x) - Hg{x)) (84) 



where both u'^{x) and TLg{x) are hermitian matrices while n^g{x) is a scalar. Furthermore, ■Ug(a;) = (^^(a;))-^, which 
implies that (x) is a real symmetric matrix. However, the conditions (73g) and (73b) require that the gluon velocity 
matrix has to be proportional to the unit matrix, exactly as that of quarks and antiquarks. 

The condition (73a) or (73c) implies that the product GGj must be proportional to the unit matrix. Therefore, the 
gluon color chemical potential must obey ll^ix) = —]lg{x). Consequently, it contains only antisymmetric components 
and it can be uniquely expressed as Jlgix) = /i"(x)T". 

Next, we analyze the conditions for cancellation of the quark-gluon collision term i.e. Eqs. (70a,70b). For SU(2) 
dabc = 0, and then it is easy to check that Eq. (70a) imposes 

T{x)=Tg{x), u^'{x)=u''g{x) . (85) 

Thus, the temperature, as well as the hydrodynamic velocity, are the same for the quark- antiquark and gluon com- 
ponents of the plasma. 

Eq. (70b) is of more complicated structure. Since it is fulfilled if = 0, let us evaluate R''. Taking into account 
that for SU(2) 

TV[tV*Q(p, x)] = '-r'^ q,{p, x) + ^S"' qo{p, x) , (86) 

and using the relation (A9) given in the Appendix, we express i?^ as 

= -lqa{p,x) TV[{T°,T^}G(pi,a;)] + ^qc{p,x) go{pi,x) - ^90(^,0;) gdPux) . (87) 

And now we refer to the diagonal gauge where the quark chemical potential is of the form 'jl{x) = ij,3{x)t^. Requiring 
= = implies gi = 92 = 0, which, in turn, demands that the respective components of the gluon chemical 
potential vanish i.e. Hg = Hg = 0. Demanding i?'^ = is only fulfilled if 

M3(a;) = fJ.l{x) . (88) 

Thus, Eq. (70b) is satisfied in arbitrary gauge if the relation (44) holds. 

The dominant processes that have been considered till now do not introduce any relation between the quark and 
antiquark baryon chemical potentials and they do not constrain the scalar gluon potential It is not surprising as 
these processes do not change the number of quarks, antiquarks or gluons. To get the relation between /i;,, fJ-b and 
/Lig, the subdominant process of quark- antiquark creation or annihilation has to be taken into account. Let us analyze 
this process. The color structure of Eqs. (76a, 76b) is rather complex. However, one checks that these equations are 
solved by the local equilibrium function (42) in the diagonal gauge (47,48). In particular, one finds that 

Tr[QQi] + 4Tr[Q] Tr[Qi] = e-^^''<P+P^^-^"'+^^''^ (lO + 4 e/^^^ + 4 e-'^^'^^ , (89) 
2Tr[G''^G"i] - 2Tr[G'G"i] + 2Tr[G"] Tr[G'i] = e-'5("-(p'+Pi)-2''p) (iQ + 4 e'^''' + 4 e"''^^) . (90) 

Thus, Eq. (76a) demands 

/Xfe + Mb = 2fi°g . (91) 

While the checking is rather simple for Eq. (76a), it is much more difficult for Eq. (76b). To reach the goal we have 
expressed the (anti-)quark and gluon distribution functions through the projections (17) and (23), respectively, and 
we have used the formula (24). Then, one finds that Eq. (76b) is satisfied if the relation (91) holds. 
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To get the chemical potentials as in the local equilibrium function (42) the binary processes are insufRcient. One 
easily observes that the equilibrium with respect to the process gg ^ ggg implies /i^ = 0. And then, Eq. (91) provides 

In summary, the requirement of equilibrium with respect to the dominant binary processes provides the local 
equilibrium state with the color structure as that in (42) which comes from the collisional invariants. The (scalar) 
chemical potentials of quarks, antiquarks and gluons are, however, independent from each other. To get the relations 
/ig = and jib = —jib, the multigluon processes and antiquark-quark annihilation into gluons must be taken into 
account. This means that the local chemical equilibrium is reached at longer time scale than the color equilibrium. 



D. Local Equilibrium Solution for the SU(iVc) plasma 



We find here the local equilibrium solutions for the SU(A^c) plasma. The quark-quark and antiquark-quark scattering 
processes are treated as for the SU(2) case. The solutions of Eqs. (66,68) read 



Q{p, x) = exp - l3{x) {u^{x)p'' - Hb{x) - ii{x)) 
Q(p, x) = exp - I3{x) {u^{x)p'' + p,b{x) + /u(a;)) 



(92) 
(93) 



The conditions for cancellation of the collision terms discussed in Sec. IV A, which involve gluons, are much more 

complicated than those for SU(2). Here, we will treat them perturbatively only. The requirement of vanishing of the 
collision term representing gluon-gluon scattering is expressed by Eq. (71). We first note that a distribution function 
proportional to the identity matrix, which is of the form exp[ — f3g(^UgPi, — /U^)] satifies this equation. We now look 
for more general solutions written as 



G{p,x) = exp - Pg{x){u''g{x)p^ - fjPg{x)) F[a{x)] , 



(94) 



where we have factored out the U(l) part of the distribution function; F is an arbitrary function of a{x) = Pg{x)Jlg{x) 
with Jlg{x) being any hermitian {N^ — 1) x {^c ~ 1) matrix. Prom Eq. (71) one deduces that F should obey the 
quadratic equation 



(T''F[a]T^ - ^ {T^'T", F[a\}^Tr[F[a\T^T''] = 



(95) 



which is trivially satisfied by the unit matrix. We now assume that F allows for an infinitesimal expansion in a 
around the identity. Then, 



F\a] 



l + a + 



and Eq. (95) imposes 



T"[5,T"] -I- -T''Ti[T'=a] = . 



(96) 



(97) 



If a is proportional to the unit matrix, the equation is obviously satisfied. However, we exclude this possibility since 
a scalar chemical potential was already included in the U(l) part of Eq. (94). A different solution of the equation is 
given by a = aoT"-- With the last option, we solve Eq. (95) to second order in 5, and find 



F[a\ = l + a+^ + 



(98) 



In principle, one can solve the equation iteratively order by order in a but the procedure becomes more and more 
difficult with every order. We will not follow it but the above results suggests that the general solution is of the form 



G(p, x) = exp - I3g{x){u''g{x)p^ - iJ,°{x) - Hg{x)) 



(99) 



as it should reduce to the Nc = 2 solution (84) with the scalar hydrodynamic velocity. 

We now look for the quark-gluon scattering, and solve Eq. (69) perturbatively for small color chemical potentials 
of quarks and gluons. In 0-th order, Eq. (69) imposes 
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T{x) = T,{x) , u'^ix) = 



(100) 



In the first order in the color chemical potentials, we find that these should obey 

(r-flixy - i {tV,/I(x)} )Tr[r^r«] + [r", T'']Tr[/I,(x)T''n = , (101) 
which is only satisfied if 

Jlg{x) = 2T»TV[/I(a;)r»] . (102) 

One could go to higher orders in the expansion but the procedure becomes very tedious. 

In the same way, one can treat the remaining processes such as the quark-antiquark annihilation. They lead to the 
same constraints as those for the SU(2) plasma expressed, in particular, by Eq. (91). 

The perturbative treatment presented here is concluded as follows. At zeroth order, the various collision processes 
allow one to fix the variables which are scalar in color space - the temperature and hydrodynamic velocity. At first 
order, every collision term imposes restrictions on the form of the matrix chemical potentials. Solving the conditions to 
all orders should simply provide the solutions, which for classical statistics are exponential functions of color chemical 
potentials. 



V. CHROMOHYDRODYNAMICS 



The form of the local equilibrium distribution function determines the character of hydrodynamics obeyed by the 
system. Here, we are going to discuss the hydrodynamic equations corresponding to the local equilibrium state found 
in the previous sections. As we have shown in Sec. IV C, the dominant processes, which are responsible for establishing 
the colorful equilibrium, do not equilibrate the system with respect to the scalar chemical potentials. The relations 
/U° = and fj,b = —fib are achieved at longer time scales. Since we are mostly interested here in the role of color 
charges in the hydrodynamic evolution, we neglect complications caused by the lack of chemical equilibrium and we 
use the distribution functions (42) where the relations /x^ = and fit = —fib are built in. 

The equations of hydrodynamics are provided by the macroscopic conservation laws of the baryon charge (32), 
energy-momentum (33) and of the color charge (34). Substituting the local equilibrium distribution functions (42) 
into Eqs. (7,8,5), one gets the baryon current, the energy- momentum tensor and the color current which enter the 
equations of ideal hydrodynamics where dissipative effects are neglected. These quantities read 



bf'ix) =6(x) u^ix) , 

t'^'ix) = [eix) +p{x)] tt^(x) u''{x)-pix) g'"' , 
j^'ix) = p{x) u'^ix) , 



(103a) 
(103b) 
(103c) 



where b, e and p are the densities of, respectively, the baryon charge, energy and color, while p denotes the pressure. 
In contrast to b, £ and p which are color scalars, the color density p is a, Nc x Nc matrix. All these thermodynamical 
quantities are given as 



2NfT^ 
37r2 



£ = 3p = 



QT* 



P = -9 — 



Nfle^^" Tr[e^l^] + e-f^^" Tr[e-'^'^]) + Tr[e'^'^''] 
/ (e'^'"' (e^'^ - Tr[e'^'^]) - e''^^"' {e'^'' - Tr[e-'^'^])) + 2T"Tr[T"e'^'^«] 



(104) 
(105) 
(106) 



Now, we consider Eq. (33) representing the energy-momentum conservation. It is well known [29] that projecting 
the continuity equation of the energy-momentum tensor on the hydrodynamic velocity, one gets the condition of the 
entropy conservation during the fluid motion. Let us see how it works here. Multiplying Eq. (33) by u'', we get 

= , (107) 

because m^m'* = 1 and u^u,^F^" = 0. The latter equality holds due to the antisymmetry of F^". Eq. (107) gives 
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(108) 



which can be rewritten as 

Ta^(sw'') + ^lb^^,{}m^') + ^v[]l^^,{fm^')] = , (109) 
by means of the thermodynamic relations 

de = Td.s + fi,4b + dp] , (110) 
e+p = Ts + fibb + Tr[/Ip] , (111) 

where s is the (local) entropy density in the fluid rest frame. The second term in Eq. (109) vanishes due to the 
conservation of the ideal baryon flow (103a) and the third term also vanishes as 

TiiJld^ipu^')] =. TviJlD^ipu'')] = . (112) 

The first equality holds because pt and p commute with each other, and consequently Tr[/i[A'^,p]] = 0. The last 
equality expresses the covariant conservation of the ideal color current (103c). Thus, Eq. (109) finally gives the 
entropy conservation 9^(.sm^) = 0. 

The analog of the Euler's equation is obtained from Eq. (33), projecting it onto direction perpendicular to u^. 
Equivalently, we consider the following combination of Eqs. (33,107) 

dj"'' - u'^u.d^t"^ = 2Ti[j,F''^] , (113) 

which gives 

(e + p) u"d,u^ = (5" - w''w,5")p + 2Tr[j,F'^"] . (114) 

To get a better insight of the physical meaning of Eq. (114) we write it down in the three-vector notation where 

n'' = (7c,7v), f = {cp:i), F°' = E\ F'^=eijkB\ (115) 

with 7 = (1 — v^/c^)~^/^ and E, B being the chromoelectric and chromomagnetic field, respectively. We have 
restored here the velocity of light c to facilitate taking the nonrclativistic limit of the derived hydrodynamic equation. 
Subtracting Eq. (114) for = multiplied by f'/c from Eq. (114) for p, = i, one gets 

(I + ^ = - + ^ 4> - ^^[^^ 4 v(j . E) + 1 j X B] . (116) 

which in the nonrclativistic domain (ti^ -C c^) reads 

(£ + rt + vV)v = - Vp - 2Tr[pE + ^ j x B] . (117) 

We note that the nonrclativistic limit, which is taken for the sake of comparison with the analogous equation of the 
electron-ion plasma [1], is only applied to the hydrodynamic velocity. The motion of the fluid's constituents remains 
relativistic. 

Although the quark-gluon plasma is composed of partons of several colors, the hydrodynamic equation (114) 
describes a single fluid. This happens because there is a unique hydrodynamic velocity in the local equilibrium state 
(42). Various color components, which enter the energy- momentum tensor, do not neutralize each other but they 
are 'glued' together in the course of evolution. Such a single fluid chromohydrodynamics was briefly considered long 
ago [6] within kinetic theory. An equation very similar to Eq. (114) has been recently derived [10] directly from a 
postulated Lagrange density. The color current, which enters the Euler's equation discussed in [10], is of the form 
QJ'^ where Q is the color charge and J'^ is the conserved Abelian current. As seen in Eqs. (103), J''' can be identifled 
with the baryon flow when we deal with a system of quarks only. In a multi-component plasma, however, such an 
identiflcations is not possible because vanishing of the baryon current does not imply vanishing of the color current. 
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VI. DISCUSSION AND SUMMARY 



Local equilibrium is only a transient state of a non-equilibrium system in its course towards global equilibrium. Thus, 

the question arises how fast such a state is achieved, and for how long it survives. We denote the two characteristic 
times of interest as tq and ri . As we have shown in Sec. IV, the dominant processes, those with the soft gluon exchange 
in t or u channel, are responsible for establishing the colorful equilibrium. Since the electric forces are screened at 
momentum transfers smaller than the Debye mass (mjj) the largest contribution to these processes comes from the 
small angle scatterings due to the magnetic forces which are effectively screened at momentum transfers below nio- 
We identify To with the relaxation time related to such interactions. Then, according to the estimate [19] found for 
the quark-gluon plasma in global (colorless) equilibrium where m£) ~ gT, we have 

-^g^T\n{l/g). (118) 

To 

We note, however, that the relaxation time in the colorful plasma can significantly differ from (118) due to the 
interaction with the background chromodynamic field generated by the color current (45) . 

And for how long does the colorful equilibrium exist? The answer crucially depends on the process which is 
responsible for the plasma neutralization. We have explicitly shown that the dominant processes comply with the 
finite color chemical potentials. We have also checked that equilibration with respect to the process qq <-> gg leaves 
the system colorhil. We expect that the collisions, even those beyond binary approximation, do not demand vanishing 
of the color chemical potentials. The point is that in every collision process, which changes the particle momenta but 
not their 'macroscopic' positions, the color current is conserved. Therefore, the collisions do not alter the local color 
charge. 

The plasma is presumably neutralized due to the collective phenomena: dissipative color currents and damp plasma 
waves both caused by uncompensated color charges. Then, the characteristic time of the system neutralization is 
controlled by the color conductivity which is again related to the estimate (118) [19]. Thus, the two times of interest 
To and Ti are of the same order, and a much more careful analysis is needed to establish the domain of applicability 
of the local equilibrium solution found here. Such an analysis should take into account not only the interaction with 
the background fields present in the colorful equilibrium but the initial non-equilibrium configuration should be also 
specified. 

At the end let us summarize the most important results of this study. The local equilibrium state dictated by 
the coUisional invariants, which follow from the energy-momentum, baryon number and color charge conservation, is 
colorful i.e. there is a non- vanishing color current in the system. The baryon chemical potentials of quarks and of 
antiquarks and the scalar (colorless) chemical potential of gluons are constrained as in a global equilibrium: fib = —^b 
and /ig = 0. The local equilibrium configuration resulting from the cancellation of collision terms, which represent 
the most probable binary parton interactions, is also colorful with the same color structure. The colorless chemical 
potentials, however, are unconstrained. The global equilibrium relations among them emerge when the subdominant 
processes are taken into account. It is conjectured that not only binary but even multi-parton collisions comply with 
the finite color chemical potentials, thus suggesting that the color neutralization of the plasma occurs not due to the 
collisions but due to dissipative collective phenomena. Proper identification of these processes and their quantitative 
description will be very important for understanding of the whole equilibration scenario of the quark-gluon plasma. 
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APPENDIX A: EVALUATION OF TRACES 

We collect here some useful formulas of the traces computed both in the fundamental and adjoint representation. 
Due to the identity 

T°-.Ti = 1 - , (Al) 

ij M 2 2Nc 
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we have the relations of the traces in the fundamental representation 



TV[rM t"B] = -J^Jr[AB] + i Tr[A] Tr[B] , (A2) 
Tr[rM] Tr[r^B] = \ Tv[AB\ - ^ Tv[A] Tr[B] . (A3) 



Furthermore, from Eq. (Al) one can deduce 



rVV" = -—t" (A4) 



2Nc 2N, 
Taking into account 

11 i 

TaU = '^Jf^ab + -^dabcTc + -^fabcTc (A5) 

one evaluates traces of products of generators in the fundamental representation. In particular, one finds 

TtIt^t^t"] = J (d''^'^ + ir^'^^ , (A6) 

_|_ _^ (^^^^ jcdr fjf^^''^ jbdr _|_ ^o,dr j^bcr^ (A7) 

For Nc = 2 one also has 

Tr[T"T''r'=r''r^] = ^ (^^"■'^fcd ^ ^cd^abe _^ ^bd ^aec _^ ^bc jade^ (-^g-. 

The identity analogous to (Al) for the adjoint representation of the SU(2) group is 

TbcTde = 5^^^''^ - S^'^^"^ , (A9) 

and we have the following relations 

Tr[T"-A T^B] = Ti[A] Tr[B] - Tr[ylB^] , (AlO) 
Tr[T"yl] Tr[T"B] = Tv[AB] - Tr[AB^] . (All) 

Using the identity (A9) one also finds 

Tr[T"T^T^] = if'"' , (A12a) 

TrlT'T^-T'T'] = S"'^'S"^ + S°'^S'" , (A12b) 

rp^^arpbrpcrpdrpe'^ ^ad jecb ^cd jeab ^ab jecd (A12c) 

In the adjoint representation of SU(3), we have the identity 

TbcTde = {S'^S"" - 5''' 5^^) - (dMfdcef ~ dbefdcdf) , (A13) 

which, in particular, allows one to compute the totally symmetric trace of 4 generators as 

iTr[{T„, Tfe}{Te, Td}] = 1 {26''%^ + S'^'Sm + S'^Xc) + Jrfab.rfcd. • (A14) 
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